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Aquifer dilation from shock waves produced by the 8.4 magnitude Alaskan earthquake of 1964 led to water 
and sediment ejection from the ground up to 400 km away from the earthquake's epicenter [1]. Groundwater 
disturbances were observed as far away as Perry, Florida ( ~5500 km), where well water fluctuations with an 
amplitude of as much as 2.3 m were reported [2]. The martian cratering record provides evidence that the planet 
has experienced numerous seismic events of a similar, and often much greater, magnitude. Given this fact, and the 
photogeologic evidence for abundant water in the early crust, we have investigated the response of a basalt aquifer 
to the propagation of compressional waves (P-waves) produced by impacts in the 33 - 1000 km diameter size range. 
The resulting one-dimensional changes in effective stress and pore pressure were calculated — as a function of both 
distance and time - based on the following assumptions: (i) that all of the seismic energy radiated by an impact is 
transmitted as a single compressional wave, (ii) that both the host rock and groundwater are compressible, and (iii) 
that there is no net flow between the water-filled pores. 

After Croft [3], an impact producing a crater of final diameter D* t has a corresponding maximum 
transient diameter D fc , given by 

D tc= J D c ° 15+a04 D/ ) - 85+004 (1) 

where D c is the transition diameter between simple and complex crater morphology, which occurs on Mars at a 
crater diameter of ~ 6 km. According to Grieve and Cintala [4], the relation between the kinetic energy of the 
impactor and D ic is given by 
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where U is the velocity of the impactor, taken as 10 km s' 1 , and g is the acceleration of gravity. However, upon 
impact, only a small fraction of Eg is actually converted into seismic energy, E s . After Schultz and Gault [5,6] 

E s =kE k (3) 

where the seismic efficiency factor, k, is taken to be equal to 10" 4 . 


In this analysis, the propagation of the resultant seismic wave is represented by a sine wave with a period 
t = 4*f 0 , where t c is the time of formation which, after Schmidt and Housen [7], is given by 

t 0 = 0.62 £ fc °-13 C/-0.04 £-°. 61 (4) 


This quantity represents the elapsed time from impact until the wavefront reaches the maximum radius of the 
transient cavity, r Q (= D tc / 2) . The maximum stress of the P-wave occurs at r Q and is equal to 
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where p is the density of the martian crust and c is the velocity of the P-wave given b> 




( 6 ) 


where E is the Young's modulus(= X + 2n, whereX andp are the Lame parameters) [5]. After Schultz and Gault 
[6], the maximum stress associated with the outward propagating shock wave is given by 


a 


a ^ n 


(7) 


The displacement of a particle at this distance, is assumed to oscillate according to the equation 


u = A q sin( k x /) (8) 

where o> is the angular frequency (= 2 jt/t), k is the wave number (= oj /c), and where, after Jaeger and Cook [8], 
A Q is the maximum amplitude given by 
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The change in strain, dc, is then calculated from the derivative of eq. (8), such that 


( 9 ) 


ds = A Q k cos(/c x - a > t) 
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After Pande et al. [9], the change in pore pressure, dp, is related to the change in strain, dc, and change in 
effective stress, da', by 

dp = K f (de - (11) 

where 


K f 
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and where n is the porosity of the material, K s is the bulk modulus of the solid matrix, and K w is the bulk modulus 
of the water saturating the pores. 


The change in strain is related to the changes in total and effective stress (do and do ' respectively) through 
the modulus matrix D* and drained modulus matrix/) [9], where 


<to =D*de 
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and where D* and D are related to each other by 
D*-D + Kf- KtD 
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(15) 


For the one-dimension case, D* reduces to E and D can be solved from eq. (15). Substituting these values 
in eq. (11), the final equation for the change in pore pressure is found to be 


To place the calculated pore pressure changes in perspective, note that during the Alaskan earthquake of 
1964, changes on the order of 1.7 bars were produced in silt and clay sediments at distances up to 400 km away 
from the earthquake's epicenter (the maximum distance at which water and sediment ejections were observed). 

On Mars, an impact of equivalent seismic energy (Dj = 33 km [10]) will produce this same change in pressure in 
basalt at a distance of ~ 100 km, while impacts with diameters > 500 km will produce pressure changes in excess of 
1.7 bars on a global scale. Impacts with final diameters > 1000 km are capable of generating pore pressures in 
excess of 10 bars out to distances of over 2000 km. Given a more realistic representation of seismic wave 
propagation through a planetary body, these pore pressures changes are likely to be amplified enormously as the 
seismic waves converge at the antipode. These results suggest that seismic disturbances produced by large impacts 
may have played a role in triggering the martian outflow channels. For example, under conditions where the local 
hydraulic head in a confined aquifer is already near lithostatic levels, the excess pressure generated by a major 
impact could conceivably disrupt a several km -thick layer of frozen ground over global distances — permitting the 
catastrophic discharge of the accumulated reservoir of groundwater until the local hydraulic head declined to the 
level of the surrounding topography. Efforts to extend this analysis to include a more realistic representation of 
seismic wave production and propagation, and to consider the effects of flow in a heterogeneous aquifer possessing 
both intergranular and fracture porosity, are currently underway. 

References: [1] Waller, R.M. (1968). in: The Great Alaska Earthquake of 1964 . National Academy of Sciences, 97-116. [2J Cooper, H.H., J.D. 
Bredehocft, I.S. Papadopulos, and R.R. Bennet (1968). in: The Great Alaska Earthquake of 1964 . National Academy of Sciences, 122-132. [3] 
Croft, S. K., (1985)./. Geophys. Res. , C828-C842. [4] Grieve, FLA.F., and M.J. Cintala. (1992 ).Meteohtics, 27 , in press . [5] Schultz, P.H., ami 
D.E. Gault. (1975a). The Moon, 12, 159-177. (6] Schultz, P.H., and D.E. Gault. (1975b). Proc. Umar Sci. Con f. 6th. , 2845-2862. [7] Schmidt, 
R.M., K.R. Housen (1987). Int. J. Impact Engng. 5, 543-560. [8] Jaeger, C., and G.W, Cook (1976). Fundamentals of Rocks Mechanics , 
Hapman and Hall., 352-254. [9] Pande, G.N., G. Bear, J.R. Williams. (1990). Numerical methods in rock mechanics . John Wiley & Sons Ltd., 
143-145. [10] Clifford, S. (1992). /. Geophys. Res., in press. 


